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Abstract 

A gap in phase space is opened up by a binary supermassive black hole as it ejects stars in a galactic nucleus. 
This gap must be refilled before the single black hole that subsequently forms can disrupt or accrete stars. We 
compute loss cone refilling rates for a sample of elliptical galaxies as a function of the mass ratio of the binary 
that preceded the current, single black hole. Refilling times are of order 10 10 yr or longer in bright elliptical 
galaxies. Tidal flaring rates in these galaxies might be much lower than predicted using steady-state models. 
Subject headings: stellar dynamics, galaxies: nuclei, black holes 
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1. INTRODUCTION 

The loss cone of a black hole (BH) at the center of a galaxy 
is defined as the set of orbits that intersect the BH, or that pass 
within some distance of its center. For instance, the tidal dis- 
ruption loss cone consists of orbits with pericenters below r t , 
the radius at which tidal forces from the BH would disrupt a 
star. Stars on such orbits are removed in a single orbital period 
or less, and subsequent feeding of stars to the BH requires a 
re-population of the loss cone, which is typically assumed to 
be driven by gravitational encounter s between stars. 

Classical loss-co ne the ory (iBahcall & Woil 119761 

Light man & Shapiro! 1 19771 ICohn & Kulsrudl 119781) was 
directed toward understanding the observable consequences 
of massive BHs at the centers of globular clusters. Globular 
clusters are many relaxation times old, and this assumption 
was built into the theory, by requiring the stellar phase space 
density near the BH to have reached an approximate steady 
state under the influence of gravitational encounters. In a 
collisionally-relaxed cluster, the density of stars near the BH 
has the Bahcall-Wolf (1976) "zero-flux" form, p oc r~ 7 / 4 , and 
the dependence of / on J, the stellar angular momentum per 
unit mass, near the loss-cone boundary is described by the 
Cohn-Kulsrud (1978) boundary conditions. The feeding rate 
is determined by the gradients of / with respect to J at the 
loss cone boundary and by the normalization of / at values of 
J far from the loss cone, i.e. by the stellar density. 

Galactic nuclei differ from globular clusters in that the re- 
laxation time is of ten very long, typically in excess of 10 10 yr 
<Fahereta1.IM997i) . One consequence is that the stellar den- 
sity profile near the BH need not have the Bahcall-Wolf r~ 7 / 4 
form. But the fact that galactic nuclei are not collisionally re- 
laxed also has implications for the more detailed form of the 
phase space density near the loss cone boundary, and hence 
for the feeding rate. In one widel y discussed model for the 
forma tion of galactic nuclei (Begelman. Blandf ord. & Reesl 
1980), a binary BH forms following the merger of two galax- 
ies. The binary ejects stars on orbits such that J < Jbi n k, 
\J2GM\2ah, with Mn the mass of the binary and a/, G/j/4a 2 
the binary's semi-major axis at the time of its formation; 

= M\MilM\i is the reduced mass and a is the stellar veloc- 
ity dispersion. Since a„ ^> r,, there will be a gap in angular 
momentum space around the single BH, of mass M. = M12, 



that subsequently forms, corresponding to stars with J <; Jj,i n 
that were ejected by the binary. Before the single BH can be- 
gin to consume stars at the steady-state rate, this gap needs 
to be re-filled. This argument suggests that the feeding rates 
of BHs in galaxies that formed via mergers might be much 
lower than predicted by application of the steady-state the- 
ory. Here we solve the time-dependent equations describing 
the evolution of the stellar phase-space density around a BH 
in a galactic nucleus and compute the time required for the 
steady-state feeding rate to be reached. 

2. approximate timescales 

Diffusion of stars into a single BH is dominated by scatter - 
ing onto low-angular- momentum orbits ( Frank & Reesll976l) . 
The scattering time is 



r 9 (r)«9(r) 2 r r (r) 



(1) 



where 8(r) is the angle within which a star's velocity vector 
must lie if it is in the loss cone, and T r is the relaxation time, 



v ; 7tG 2 w^p(r)logA' 



(2) 



with ra* the stellar mass and log A » 15 the Coulomb loga- 
rithm. The square-root dependence of 8 on Tq reflects the fact 
that entry into the loss cone is a diffusive process. If the loss 
cone was initially emptied of all stars with pericenters below 
some radius ro, the time to diffusively refill this angular mo- 
mentum gap is 

(3) 



Equation[3]assumes that r, <C ro S5 fh with r/, = GM\ ? / G 2 the 
influe nce radius of the single (coalesced) BH ( Frank & Rees 
[I976T) . Until the angular momentum gap is filled, the rate of 
supply of stars to the BH's capture sphere will be much less 
than its steady-state value. 
In the binary coalescence model, most stars with pericenters 
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r p <Ka h =K^ 



(4) 



will have been ejected prior t o coalescence by the grav- 
itation al slingshot mechanism (Saslaw, Valtonen & Aarseth 
with K a constant of order unity. This ejection will 
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occur even if the binary's main source of angular momentum 
loss is torques from gas clouds, as long as \a/d\ exceeds stel- 
lar orbital periods near the BHs. Setting ro w r p and r w r% 
(since most of the scattering into the single BH takes place 
from radii near r&), we find 

Tga ^l{TTW Trirh) (5) 

withq=M 2 /Mi < 1. Foro = 1, T gap w T r /16, and for g< 1, 
Tgap ~ (?/4)r r . These times are short compared with T r but 
still in excess of 10 10 yr for the luminous galaxies that are 
most likely to have experienced mergers. 

The principal uncertainties in this estimate of the refilling 
time are the value of K, and the fact that stars are scattered into 
the loss cone from a range of distances with different values 
of T r . We address these issues by more careful calculations in 
the following sections. 

3. CREATION OF A PHASE-SPACE GAP 

The gap is created by the binary BH when its separation 
reaches ~ a/,. The definition given above for is difficult 
to apply to galaxies in which a is a function of radius. We 
therefore adopt an alternative definition in terms of r/,: 

with r/, defined as the radius in the unperturbed galaxy con- 
taining a mass in stars equal to twice M\2- These defi- 
nitions are equivalent to the standard ones r/, = GMn/o 2 , 
a/, = G/j/4o 2 in a singular-isothermal-sphere nucleus but can 
also be applied to nuclei in which a is a function of radius. 

Once it forms, the binary shrinks as it ejects stars that pass 
within a distance^ a(t) of the binary's center of mass. The 
coupled evolution of the stellar fluid and the binary is com- 
plicated due to processes like re-ejection, the repeated inter- 
action of (non-escaping) stars with the binary, and by the fact 
that th e binary's cross-section is changing with time (Merritt 
2004). 

To model this process, and to compute the size of the gap 
in terms of a/, as we have defined it, we carried out a set of 
Monte-Carlo simulations. Scattering experiments were first 
carried out for an isolated, circular-orbit binary as a function 
of the parameters (R p ,V p ), the distance and velocity at closest 
approach to the binary that a test particle would have if the bi- 
nary were replaced by a point of mass M\2 at the binary's cen- 
ter of mass. The changes in energy and angular momentum, 
AE and AJ, experienced during the encounter were recorded 
and the calculation repeated for a large number (~ 10 6 ) of 
different (R p ,V p ) values. 

In the second step, a set of Monte-Carlo positions and 
velocities were generated for a set of stars drawn from an 
isotropic distribution function describing the stars in the un- 
perturbed galaxy. The gravitational potential was fixed as 
<E>(r) = 4>*(r) — GM\ijr with <E>*(r) the contribution to the 
potential from the stars. The stellar orbits were integrated 
forward in time, and stars were given kicks in energy and 
angular momentum whenever they passed within a distance 
r C rit(t) = 3a(f) of the galaxy center. The kicks were drawn 
randomly from the stored values of AE and AJ ; the star's or- 
bit in the full potential was related to the Keplerian orbit in 
the scattering experiments by equating the pericenter distance 
r p and velocity at pericenter v p of the orbit about the galaxy's 
center with R p and V p . Immediately after a kick, a star was 
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FIG. 1 . — Lindblad diagram of stars at the final time step in the Monte- 
Carlo experiments, with q = 1 (a) and q = 1/8 (b). Additional parameters are 
given in the text. Larger dots are stars that are still on orbits that intersect the 
binary. Lines show j gap for K = 0.5 (blue), 1 (black) and 2 (red). 



placed on an outgoing orbit on the sphere r = r cr u with a ve- 
locity computed using its new values of E and J. At the same 
time, the binary's energy was reduced by m+AE and its semi- 
major axis correspondingly changed. Possible changes in the 
eccentricity of the binary's orbit were ignored. Stars with ini- 
tial apocenters within a(Q) were flagged and removed; in a 
real galaxy, such stars would have been ejected during the bi- 
nar y's formation. T he stellar potential was taken to be that 
of a lDehnenl dl9 93) model; initial positions and velocities of 
the stars were generated from the unique isotropic distribution 
function f(E) that reproduces the Dehnen mass distribution in 
the combined po tential of the stars and the central point mass 
llTremaine et alJlT994l) . The parameters of the Monte-Carlo 
integrations were q, M\2/M ga [, a(0), and y, the central loga- 
rithmic density slope of the Dehnen model. 

Figure[Oillustrates the gap created in two integrations with 
y = 2 and M n /M gal = 0.001, and with q = (1, 1/8). The plots 
show the phase-space distribution at a time of ~ 4 in units 
such that G — M ga i = rjy = 1 with ro the Dehnen-model scale 
length. The larger circles are stars that are still interacting 
with the binary, i.e. stars with r p < r cr i t . Most of these stars 
will eventually be ejected but some will remain in the galaxy 
as the binary "shrinks away" beneath them. The three curves 
in FigureQare 

J 2 = 2Ka h [E - <S>{Ka h )] = J 2 gap (7) 

with K = (0.5, 1,2); J sap is the angular momentum of a star 
with pericenter Ka/,. There is a sharp drop in phase space 
density at J w J gap, K w 1. This relation was found to hold 
also for other values of y and q. The location of the gap was 
found to be only weakly dependent on the choice of a(0). For 
the integrations of Figure^ a(0) « 2a/,, and increasing a(0) 
by a factor of 10 caused J gap to increase by only ~ 50%. We 
conclude that equation Q with K = 1 is a good representation 
of the phase-space gap produced by a binary BH. 

4. REFILLING THE GAP 

In the Monte-Carlo experiments described above, the bi- 
nary decay stalls after all the stars on orbits intersecting the 
binary have been ejected. We assume that some other mecha- 
nism, e.g. torques from gas in an accretion disk, then acts to 
extract angular momentum from the binary, eventually bring- 
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ing the two components close enough together that emission 
of gravitational radiation can induce complete coalescence. 

The phase-space gap opened up by the binary will then 
gradually refill on the time scale T gap . We modelled this evo- 
lution using the time-dependent Fokker-Planck equation. The 
following simplifying assumptions were made. (1) Changes 
in stellar energy due to encounters were ignored, i.e., we as- 
sumed that Tg al , -C T r (cf. equation|5}. (2) Relevant stars are 
those with R = J 2 / J C {E) 2 -C 1, where J C (E) is the angular 
momentum of a circular orbit of energy E. Under this as- 
sumption, the low-/? form of the velocity diffusion coefficient 
may be used. (3) The diffusion in / is slow compared with 
orbital periods so that the orbit-averaged form of the Fokker- 
Planck equation can be used. We justify this assumption be- 
low. Given these assumptions, the evoluti on of the s tellar 
phase space density f(J,t;E) obeys ( Milos avlievic & Merittl 
2003J 

3/13/ .3 A .„ 
3^-4737^37; (8) 

where j = J/J C (E) = R l l 2 is a dimensionless angular mo- 
mentum variable and x = a(E)t is a dimensionless time, with 
a(E) the orbit-averaged diffusion coefficient at energy E: 

((AR) 2 ) 



a(E) 



P(E) 



dr 

— lim 



2R 



(9) 



and P(E) the period of a radial orbit of energy E. Hence x « 
t/T r (E). Equation (|8j can be integrated forward in time given 
boundary conditions at j = and 7 = 1 and initial conditions 
f(j,t — 0;E). Figure 13 shows the result if f(J,t = 0) has the 
form 

f = 0, ;<0.2; / = 1, 0.2<;<1. (10) 

These are the initial conditions, at one energy, created by a 
binary BH with j gap = 0.2 in an initially isotropic nucleus. 
The boundary conditions are 

f(j<jic,t) = 0; §y|y=i = (11) 

with ji c the scaled angular momentum of a capture orbit; in 
Figure|5] ji c = 0.02. The flux of stars into the BH is given by 



N(j,t;E)df 



a dN 



(12) 



2 3logi Ufc 

with N(E,R,t)dEdR = 4n 2 P(E)J 2 (E)f(j 2 ,t;E)dEdR the 
number of stars in the integral-space element dEdR. 

The orbit-averaged assumption breaks down sufficiently far 
from the binary, in the "pinhole" regime ( Lightman & Shapiro 
119771) . In the loss cone defined by a single BH, the breakdown 
occurs at binding energies less than E cr i t where 



a (F\ a ( £ ) P ( £ ) 
9lc{E) 



■ 1. 



(13) 



Ric(E) 

For E <, E cr it, stars can wander into and out of the loss cone in 
a single orbital period and the loss cone remains nearly full in 
spite of capture . Typically r cr u as r/,, the BH's influence radius 
JFrank & Reeslll97 6l). In the case of the angular momentum 
gap created by a binary BH, the orbit-averaged assumption 
holds initially even for r ^S> r cr j t , since stars need to diffuse 
over the full angular momentum gap before they can be dis- 
rupted. Breakdown of the orbit-averaged approximation only 
occurs much farther from the center where 



_ a(E)P(E) 
q gap (E) = — — — > 1. 

Rgap(E) 




FIG. 2. — Solution to equation |SJ with the initial and boundary conditions 
given in the text, at dimensionless times 1 = 0,1- 
where X = /i(E)t fst/T r . 



l(T 3 ,2x 10~ J ,...,4x 10 



-3 



(14) 



ICohn & Kulsrudl f l978) showed that the steady-state form of 
/ near the loss cone was / ~ ln(y'/ 70) and derived an expres- 
sion for jo = jo(E); beyond r cru , jo w 0, i.e. the loss cone is 
nearly full. We adopted the same expression for jo here and 
set / = 0, j < jo at all times. While not strictly correct in the 
case of a time-dependent loss cone, this boundary condition 
gives the correct result as t — ► °°, and the error in the inferred 
loss-cone flux at earlier times should be very small since so 
little of the galaxy is in the full-loss-cone immediately after 
the binary forms. 

We applied this prescription to compute the evolution of / 
and f in the "core" galaxies from Paper I. For each galaxy, 
the function a(E) was computed on a grid of energies and 
stored. Equation ^ was then integrated forward for the ap- 
propriate dimensionless time at each E value and the flux 
computed from equation dl 2t i. We used the same values for 
the BH mass M and for r t as in Paper I (only the "M — o" 
values for M from that paper were used.) We tried three dif- 
ferent values for the mass ratio of the binary that was assumed 
to have formed the gap: q = M2/M1 = (0.1,0.3, 1), and j gap 
was computed from equation Q for various values of K; fol- 
lowing the discussion above, we focus here on the results ob- 
tained with K = 1 . At the end of the integration, the galaxy 
was assumed to be in its currently-observed state, for which 
we know f(E), the /^-averaged distribution function. The cor- 
rectness of the integrations was tested by verifying that the 
asymptotic flux was equal to the steady-state values computed 
in Paper I. 

5. RESULTS 

Figure|5]shows the energy-dependent evolution of the loss- 
cone flux in one of the sample galaxies, NGC 4168, for K = 1 
and q — 0.1. Stars with binding energies near 4>(r/,) are the 
first to be scattered into the BH, followed by stars with en- 
ergies slightly above or below £),. The fact that very lit- 
tle of the flux comes from stars with energies very different 
from Ei, validates our approximate treatment of the loss-cone 
boundary condition in the full-loss-cone regime. The total 
flux reaches (1%, 10%, 50%, 90%) of its steady-state value in 
atime of- 4.5 x 10 9 ,9.8 x 10 9 ,1.7 x 10 10 ,9.7 x 10 10 yr. 

We defined two characteristic times associated with loss- 
cone refilling, to is the elapsed time before a single star would 
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log 10 (t/Gyr) 

FIG. 3. — Dependence of the loss-cone flux on energy and time in the 
galaxy NGC 4168 assuming a binary mass ratio of q = 0.1. Curves are la- 
belled by their binding energy; the value of the gravitational potential at r = ri, 
in this galaxy is 1.76. Thick black curve is the total flux, in units of stars per 
year. 




FIG. 4. — Two characteristic times associated with loss-cone refilling in the 
sample galaxies, to is the elapsed time before the first star is captured and fj; 2 
is the time for the loss-cone flux to reach 1/2 of its stea dy-state value. Solid 
lines are the approximate fitting function, equation 1151 . 



steady-state value. Figure |4] shows values of to and t\/2 com- 
puted for our sample of galaxies, using K = 1 andg= (0.1,1). 
to is of order 10 9 yr or greater for all galaxies and for both 
values of q, and t\ n is roughly an order of magnitude greater. 
The relation: 



'1/2 



c- 



■ C- 



M. 



10 u yr 1O 8 M (1+^) 2 1O 8 M 
fits the data in Figure Intolerably well with C = 1. 



(15) 



6. DISCUSSION 

The flux of stars into the tidal disruption loss cone of 
a BH that formed via binary coalescence can be much 
lower tha n would be predicted under the steady-state as- 
sumption dSyer & UlmeJl999HMagorrian & Tremaine! 19991 
Wang & Merritt 120041) . Time scales for loss-cone refilling in 
the nuclei of bright elliptical galaxies are of order 10 10 yr or 
longer and scale roughly linearly with BH mass for a fixed 
(assumed) mass ratio of the binary that preceded the current, 
single BH. Even if the galaxy merger occurred as long as 
~5x 10 9 yr ago, the flux of stars into the BH's tidal disrup- 
tion sphere could still be less than ~ 1% of the steady-state 
value. BH feeding rates in bright elliptical galaxies were al- 
ready expected to be lower ( < 10~ 5 stars yr -1 ; Paper I) than 
in less luminous galaxies, due primarily to their low central 
densities. Based on the arguments in this paper, stellar disrup- 
tion rates in these galaxies might be much lower still, making 
the nuclei of bright ellipticals unlikely sites for observing a 
tidal flaring event. Event rates in fainter galaxies, on the other 
hand, should be closer to the steady-state values computed in 
Paper I since these galaxies are less likely to have experienced 
a merger since the epoch at which the stellar cusp was formed, 
and since relax tion times in their dense nuclei are probably 
relatively short (lLauer etalll998l) . 

We have investigated an extreme, but physically plausible, 
model in which the region around a newly-coalesced BH was 
emptied of stars out to radii far greater than the radius of the 
tidal disruption sphere. But even in less extreme nuclear for- 
mation models, the long relaxation times in stellar nuclei im- 
ply a slow approach to a steady-state distribution of stars near 
the loss cone and hence a loss cone flux that can be very differ- 
ent than the value computed via the standard, steady-state the- 
ory. Our results highlight the need for a fully time-dependent 
theory of loss cone evolution that can be applied to systems, 
like galactic nuclei, that are much less than one relaxation 
time old. 



be scattered into the loss cone. Obviously, the value of fo 
might depend rather sensitively on the choice of initial con- 
ditions. A more robust measure of the refilling time is ti/2, 
the elapsed time before the loss-cone flux reaches 1 /2 of its 
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